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Abstract. The flow curves, viz. the curves of stationary stress under steady shearing, are obtained close 

to the glass transition in dense colloidal dispersions using asymptotic expansions in the schematic F{^ - 
model of mode coupling theory. The shear thinning of the viscosity in fluid states and the yielding of 
glassy states is discussed. At the transition between fluid and shear-molten glass, simple and generalized 
Herschel-Bulkley laws are derived with power law exponents that can be computed for different particle 
interactions from the equilibrium structure factor. 

PACS. 82.70.Dd Colloids - 83. 50. Ax Steady shear flows, viscometric flow - 83.60.Df Nonlinear viscoelas- 
ticity - 83.60.Fg Shear rate dependent viscosity - 64.70.P- Glass transitions of specific systems - 64.70.Q- 
Theory and modeling of the glass transition 



1 Introduction 

The behavior of dense colloidal dispersions under flow is 
especially interesting at high concentrations, where flow 
interferes with solidification [112] . The curves of shear stress 
versus shear rate under steady flow, also called flow curves, 
provide insights into the cooperative particle rearrange- 
ments under strong external drive, and may be affected 
by yielding, shear thinning and thickening, and flow in- 
stabilities leading to heterogeneous and/or intermittent 
flow [3I4I5I6I7I8] , While the interplay between flow and 
freezing into crystalline ordered states is of interest in dis- 
persions of rather monodisperse particles, the interplay of 
flow and arrest into metastable and amorphous solids (the 
glass transition) is of relevance in more complex disper- 
sions, especially ones consisting of polydisperse particles, 
where (frozen in) disorder frustrates crystal packing [9110] . 
Flow curves thus can provide insights into strongly driven 
fluids and glasses, because the shear rate can easily be 
made much larger then the internal relaxation time To, 
and thus flow strongly affects the structural dynamics, 
viz. the strongly cooperative particle rearrangements in- 
trinsic to the glass transition. Various phcnomcnological 
formulae have been proposed to describe stationary flow 
curves, like the Cox-Merz rule, the Herschel-Bulkley law, 
and the concept of a power-law fluid [2J. Discussion of 
these phcnomcnological laws starting from a microscopic 
description of a colloidal dispersion is aim of the present 
contribution. 

The indicated connection between the glassy dynamics 
and the nonlinear rheology of dense dispersions raises the 



question for a unified description of both phenomena; see 
[8] and references therein. A microscopic theoretical ap- 
proach for the shear-thinning of concentrated suspensions 
and the yielding of colloidal glasses was presented |11I12J , 
which builds on the mode coupling theory (MCT) of ide- 
alized glass transitions developed by Gotze and cowork- 
ers |13ll4j . This description of the glass transition in a 
quiescent colloidal dispersion explains many findings ob- 
served in especially dynamic light scattering experiments 
on model colloidal dispersions |15I16I17I18I19I20I21J . and 
also gives results for the linear response viscoelastic prop- 
erties, viz. linear (frequency dependent) shear moduli. The 
extension to strongly driven stationary states [llj , the so- 
called ITT (integration through transients) approach, and 
its consecutive generalization to arbitrary time-dependent 
states far from equilibrium [22] and to arbitrary flow ge- 
ometries [23] . has yielded a non-Newtonian constitutive 
equation applicable to concentrated dispersions in arbi- 
trary homogeneous flows, albeit under the approximation 
that hydrodynamic interactions are neglected. Within the 
theory, this approximation becomes valid close to the glass 
transition and for weak but nonlinear flows, where the 
slow structural relaxation dominates the system proper- 
ties, and where hydrodynamic interactions only affect the 
over-all time scale. The structural dynamics under flow re- 
sults from a competition between local particle hindrance 
(termed cage effect) and the compression / stretching (i.e. 
advection) of the wavelength of fluctuations induced by 
the affine particle motion with the flow. 

Close to the glass transition, where the structural re- 
laxation has become much slower than any other intrinsic 
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relaxation process, ITT makes universal predictions that 
may be captured in simple, so called schematic models, 
where a few parameters mimic the microscopic input into 
the ITT approach [24] . These so-called vertices mimic the 
equilibrium structure factor which would be required to 
solve the complete ITT equations. Some of the universal 
predictions, like the existence of power laws and their ex- 
ponents, survive the simplifications leading to schematic 
models and agree quantitatively with the microscopic ITT 
results for e.g. hard spheres. Results from one of the sim- 
plest schematic models, the F^ -model, have been com- 
pared to flow curves measured in dense dispersions of col- 
loidal hard spheres [25 10 8 2 8|29] and in simulations of 
hard spheres and supercooled binary Lenard-Jones mix- 
tures [26|5|27| . The model explained salient points of the 
flow behavior of glassy systems, like a finite dynamic yield 
stress beyond the glass point. It vanishes discontinuously 
when going below the glass transition, where a first New- 
tonian plateau appears in the viscosity, which is followed 
by strong shear thinning. 

Asymptotic expansions capture the leading universal 
results of the schematic model and allow to discuss the 
complete scenario of the nonlinear rheology as function 
of shear rate and separation to the glass transition. They 
thereby provide insights into the connections between var- 
ious rheological features as described by ITT. The phe- 
nomenological laws of flow curves can be discussed. In 
this paper we present the asymptotic solution of the F^ - 
model, thereby supplementing its use when fitting exper- 
imental data. The model is summarized in Sect. 2. The 
nonlinear stability analysis, determining the existence of 
power laws and their exponents is presented in Sect. 3. It 
is this stability equation which holds in the microscopic 
ITT equations as well, and which provides the link to cal- 
culate e.g. the exponents appearing in flow curves from 
first principles within ITT. Sect. 4 discusses the yield pro- 
cess, where the structural relaxation is driven by the ex- 
ternal shearing. Sects. 5 and 6 discuss the results for the 
flow curves, which are summarized and brought into com- 
parison with experiments in Sects. 7 and 8. 



2 Schematic model 

2.1 The transient density correlator 

The schematic F^ -model considers a single typical den- 
sity correlator <P (t) obeying a Zwanzig-Mori equation [24] : 

1 



<P(t)+<P(t) + m(t-t')$(t')dt' = 0. (1) 

r jo 

r^ 1 defines the (bare) time scale of the system. A van- 
ishing memory kernel would lead to an exponential decay 
of the correlator. The equation of motion is closed by the 
self consistent approximation 



and the normalization <P (t — > 0) = 1 — Tt. Increasing the 
particle caging is modeled by increasing the coupling ver- 
tices Vi > 0. The parameter 7 is the shear rate which 
leads to a strain dependent decay of the memory kernel. 
The parameter 7 C is a scale for the magnitude of strain 
that is required in order for the accumulated strain jt to 
matter [8]. 

For a vanishing shear rate, 7 = 0, this model has been 
suggested by Gotze [30] and describes the development 
of slow structural relaxation upon increasing the coupling 
vertices Vi. In this case the model shows two transition 
lines. For sufficiently small Vi the noncrgodicity parame- 
ter / defined by the long time limit of <P (t) is zero. By 
crossing one of the transition lines, / becomes positive: 
/ = <P(t — > 00) > 0. It is proved [31] that for a van- 
ishing shear rate <P (t) is completely monotone and / is 
defined by the maximum real solution of / + g (/ — 1) = 
with g = m{t — > 00) = v\f + i^/ 2 '■ This determines the 
transition lines. Vertices lying on a transition line and the 
corresponding nonergodicity parameters are denoted by 
vi and f c . The type A transition line with f c = is given 
by vi = 1, < x>2 < 1. By crossing this line / varies con- 
tinuously. The type B transition line, where / jumps from 
zero to a finite value f c = 1 — A, can be parametrized by 
(v^v^) = ((2A — 1) , 1) /A 2 with the exponent parameter 
1/2 < A < 1. At (wf,t>§) = (1,1) the type A transition 
line merges continuously into the type B transition line. 

In the following we only consider vertices Vi = vf + 5vi 
close to the type B transition line. Then the dimension- 
less separation parameter measuring the distance from the 
transition point is denned by 



l-fc ' 



(3) 



With these definitions the model describes the dynamics 
of an ideal glass forming system. A negative value for e 
defines a liquid like state. For e —> 0~ an increasingly 
stretched plateau close to f c develops and the time scale 
for the final decay to zero depends sensitively on e. A posi- 
tive separation parameter leads to a glassy state where the 
noncrgodicity parameter is determined to leading order by 



/(e^0+, 7 = 0)=/ c + (l-/ c ) 2 



1-A 



(4) 



■(*) 



l'!<P (t) + V 2 <P 2 (t) 

1 + {ithef 



(2) 



A detailed discussion of the quiescent Fi2-model can be 
found in [T5] . 

For non-vanishing shear rates, 7 7^ 0, because of the 
vanishing long time limit of the memory kernel, also / 
vanishes for all separation parameters. An arbitrary small 
but finite shear rate melts the glass and restores ergod- 
icity. In the liquid state (e < 0), small shear rates have 
only little effect on the dynamics. But if the time scale for 
the shear induced decorrelation becomes faster than the 
time scale for the structural relaxation, then the long time 
dynamics is shear dominated and shows the same scaling 
behavior as for the glassy state (e > 0). 
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2.2 Flow curves 

In the framework of the schematic model the steady state 
shear stress a = 777, where 77 denotes the viscosity, follows 
via integrating up the generalized shear modulus [10] : 



7 / 9 (t, 7) dt, 



g{t^) = v a <P 2 {t)+T lao 8(t + {r) 



(5) 



(0) 



The small negative quantity 0~ in Eq. j6]) ensures that the 
integral in Eq. ([5]) includes the singularity of the 8- function 
completely. The quantity v a is assumed to be a constant 
and the parameter 7700 models viscous processes which 
require no structural relaxation [8j. The resulting linear 
moduli, the storage G'{u>) and the loss modulus G"(w), in 
the quiescent situation are discussed in Ref. [5] . We do not 
comment further on the Cox-Merz rule for which there is 
no basis in the asymptotic expansions of the F^ -model. 
We sketch some basic properties of the flow curves (a 
as a function of 7). In the liquid state, because of the 
existence of r/o = J Q 9 (i, 7 — 0) dt, the flow curves show 
linear asymptotes for small shear rates and the yield stress 
a + , the zero shear limit of the stress a + = a (7 — > 0), is 
zero. A simple Taylor expansion of g (t, 7) generalizes this 
statement: a = 7 (770 + 7700) + O (j 3 ) . The Newtonian vis- 
cosity 770 diverges at the glass transition according to the 
well known MCT power law, 770 oc |e|" 7 [T3]; see Eq. fig]) 
below. In the glassy state, the time scale for the shear 
induced decay scales with 1/ I7I which leads to a non van- 
ishing yield stress. At e — the yield stress jumps from 
zero to a finite value. For high shear rates, all flow curves 
converge to the common linear asymptote 7 {v a /2r + 77^) 
where v a j2F < ?] . Hence the model describes the transi- 
tion from a shear thinning liquid to a yielding glass [24] : 
In the liquid state, the viscosity starts at the plateau value 
Vo + r 7oo for small shear rates and decays monotonically to 
a second plateau value v a /2r + ^ for high shear rates. 
In the glassy state, the viscosity diverges in the zero shear 
limit. 



2.3 Parameters 

The contribution of the parameter 7700 is only an addi- 
tional constant to the viscosity which is not of interest 
for our work. Hence we choose 77^ = 0. In the following 
we choose for all numerical examples vi 
Vi = Vi + £ (1 - f c ) /f c . For 7700 =0we easily verify 



v'2 = 2 and 



*(e,7,7 ,r,t)=£ 



7 



,rt 



a (e, 7) 7c F) = v a j c a e, — - 



(7) 



(8) 



Hence, without loss of generality, we choose units in such 
a way that "f c = T = v a = I and that t, 7 and a are di- 
mensionlcss. For the numerical solution of the F{^ -model 




2 4 

Fig. 1. Numerically obtained transient correlators $ (t) (solid 
lines) for e = 0.01 (red, 1), e = (orange, 2), e = -0.005 
(green, 3) and e = —0.01 (blue, 4). All curves were calculated 
with 7 = 10~ 7 . The dashed lines show the corresponding nu- 
merically obtained functions f c + (1 — / c ) Q (t). 



we use the algorithm first published in [32] with a grid 
size of N = 256 and an initial step size of h — 10~ 6 , see 
also [33] for details. All numerical parameters, which arc 
relevant for our work, are tabulated in Appendix [Bl 



3 Glass stability analysis 

3.1 The equation of motion 

For small separation parameters and shear rates the cor- 
relators develop a stretched dynamics located around the 
critical plateau value f c . The discussion of the dynamics 
around this plateau, also called the /3-relaxation process 
[13] , is topic of this section. To derive the equation of mo- 
tion to leading order in the separation parameter and the 
shear rate we first rewrite the convolution integral in Eq. 
© as Jj m(t- if) & (f) dt' = f t /„ m (t - t') <P (£') dt' - 
m (t) . Eq. ([4]) motivates the expansion 



*(*) = /c+(l-/c)5(*) + 



(9) 



with g (t) = O (|e| 1/2 ). The function Q (t) is often called 
/3-correlator. By expanding the memory kernel in powers 
of 7 and neglecting the time derivative we obtain to lead- 
ing order the nonlinear stability equation [24] , 



c (t) (^) 2 + X g 2 (t) 



d 

di 



g{t-t')g(t')dt' (10) 



with c ( ^ = (1 - A) /A 2 . Eq. (fJO]) has to be completed with 
the initial condition 



g(t 



°>=<i 



where the exponent a obeys 

r 2 {i-g) 
r(i-2o) 



= A. 



(11) 



(12) 
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The time scale to enables one to match Eq. ([lljl to the 
microscopic dynamics in the correlators. It has to be de- 
termined numerically by matching to the full solution for 
<P(t) from Eq. |T]). For e = and 7 = 0, the power law 
occurring in Eq. (TrT| is a special solution of Eq. ([11])) and 
describes the critical decay of <P (t) to f c [13]. Fig.[T]shows 
representative examples. 



3.2 The two-parameter scaling law 

One of the important aspects of Q (t) are its homogene- 
ity properties [24] which follow from Eq. (fT0|) . With some 
arbitrary scale J? > it obeys the two-parameter scaling 
law 

G (t, e, 7 ) = Q a G (t = Q^,e= eQ- 2a , 7 = 7 t tf- (1+o) 
\ t a 

(13) 

with the master function obeying the initial condition 

Q (i — > 0) = (t) . Eq. ([T3|) allows to define three regions. 

1 
By choosing tt = J7io| 1+a we obtain 7=1 and e = e/sj 

with the natural scale 




6 8 

log 10 (t) 

Fig. 2. An overview of the properties of Q (i) (solid lines) 
for the same values for e and 7 (1-4) as in Fig. Q] The dotted 
lines show the leading asymptotes for the corresponding time 
scales: The critical decay (t/to)~ a (a), the von Schweidler law 
— {t/ro) (b), the arrest on the plateau value ye/ (1 — A) (c) 
and the shear- induced linear asymptote — t/r-y (d). The dashed 

-1.39 • 10 3 



line shows Eq. 
(e = 0). 



22} evaluated to n = 3 for at 



l7*o| 



(14) where the exponents u, v and the coefficients a, a± obey 



The liquid region is defined by e <C — e-y, where a finite 
shear rate barely distorts Q (t). In the transition region, 
|e| <C e 7 , the dynamical anomalies arc the most pro- 
nounced. The yielding glass region is defined by e 3> £ 7 . 
These definitions hold for |e| <C 1 and \ito\ -C 1 what we 
assume throughout in the following. 

The present two-parameter scaling law bears some sim- 
ilarity to the one presented by Gotze and Sjogren for the 
description of thermally activated processes in glasses [53] . 
In both cases, ideal glass states are destroyed by additional 
decay mechanisms. Yet, the ITT equations and the gener- 
alised MCT equations differ qualitatively in their descrip- 
tion of the states at e > 0. While the generalised MCT 
there predicts finite Newtonian viscosities, ITT predicts 
yield-stresses. Also both mechanisms can appear simulta- 
neously and then need to be combined [8]. The similarity 
between both two-parameter scaling laws thus underlines 
the universality of the glass stability analysis, which is de- 
termined by quite fundamental principles. In Eq. (jlOp . the 
shear rate can only be a relevant perturbation (at long 
times) if it appears multiplied by time itself. Symmetry 
dictates the appearance of (jt) 2 , because the sign of the 
shear rate must not matter. The aspect that shear melts 
the glass determines the negative sign of ("it) 2 . 



Jri) 



(jt) 2 = a 2 r Q , t 2u + 2a 2 ai r^t 2u+ \ 



(16) 



A and the other co- 



, _ r{l+u+nv)r(l+u+n'v) 

With I „,„' — r(l+2u+{n+n')v) 

efficients are determined by the recursion formula a n >2 = 

27T) 2^tn' — l n — n',n'Q'n—n'Q"n,'' 

Along the lines of the work [34] , we have analyzed all 
possible combinations that solve Eq. ([TBI) in detail, see 
[3"3"j . Each defines an (asymptotic) series expansion de- 
scribing the dynamics in a special regime on a special time 
scale. There are four solutions for the liquid region, two for 
the transition region and three solutions describing yield- 
ing glassy dynamics. Here we only present a qualitative 
overview containing the leading asymptotes. In the liquid 
region, for short times Q (t) follows (t/t )~ a and merges 
into a second power law — (t/ro) for intermediate times 
with the von Schweidler exponent b obeying 



r 2 (i + b) 
r(i + 2b) 



= A. 



The two-parameter scaling law leads to 



to = T t lei 



(17) 



(18) 



3.3 The power series ansatz 

Eq. (|10p can be solved by using a generalized power series 
ansatz 



g(t) = od u ll + J2a n t v A 



(15) 



7 



2ab 



(19) 



The parameter fo has to be determined numerically. With- 
out shear, the divergence of t from Eq. p8[) is the origin 
of the divergence of the Newtonian viscosity 



VO = G^TQTrj 



(20) 
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where G^ denotes the shear modulus (clastic constant) 
of the glass [8], and t v is a constant. For long times Q (t) 
merges into the linear asymptote —t/r-y with 



(21) 



In the transition region, after following (t/to) a the func- 
tion Q (t) merges into the long time asymptote — t/r-y. In 
the yielding glass region, Q (t) follows (i/i )~ a , arrests on 
the plateau value y/e/ (1 — A) for intermediate times and 
merges into the linear asymptote — i/r-y for long times. So 
we can summarize that the short- and long time asymp- 
totes are common for all e if 7 7^ is common. Fig. [2] 
shows an overview of the properties of Q it) . 



3.4 The long time dynamics 

We briefly present the special asymptotic series expansion 
which is essential to understand the asymptotic behavior 
of the flow curves. For long times, e can be neglected on 
the left side of Eq. (fT6|) . We solve the resulting equation 
by choosing u = 1 and Jq.i = 0. The result can be written 
as 



»(*)-— (i+E<(°i*-°r) 



(22) 



with a\ = 1, a,;> 2 = -53^ E£'=\ r n-n', „'<_„, a*, and 
the exponent c is defined by 



2A- 1 
A ' 



(23) 



A numerical example for e — is shown in Fig. O 

The coefficient a\ has to be determined by matching 
to shorter time scales [24]. It obeys the two-parameter 
scaling law 



"1 



(e,j) = &ifa- 2a ,jt Q-^)(f- 



(24) 



Because of this scaling law, the non-trivial factor of ax 
depends on the ratio e/e-y only: 



ai (e, 7) = (tof |t*o| 1+a Si (e/e-y, 1) ■ 



(25) 



This behavior will be crucial in the following as it enables 
us to describe the flow curves for small |e| and | 'yto | - Three 
limiting behaviors of a\ can be determined, by identifying 
the time scales at shorter times, to which the series from 
Eq. (|22|) needs to be matched. 



describing this crossover, we solve Eq. (|16p by neglecting 
e and choosing a = — (to) - and u — b. The resulting 

natural time scale for the crossover process is given by 

1 
T b = To I7T0I 1 ~ b to which the long time series has to be 

matched. For this, we choose J? such that t = t/rb, use 

Eq. (|18p and assume that in the limit e/e-y —>■ — 00 the 

coefficient di converges toward some constant. We obtain 



01 (e/e 



-00) = Oi (io) c |e| 



-><■- 



It*oI 



(26) 



where a x is an appropriately chosen constant. 



3.4.2 Critical point 

At the glass transition point, Eq. (|25[) leads directly to 



ai 



(0, 7 ) = a?(to) c |7to| 



(27) 



with some constant a\ 



A matching to the natural time 
scale, where the crossover from the critical short time dy- 
namics to the shear-induced long time dynamics occurs, 
leads to the same result. 



3.4.3 Yielding glass 

In the yielding glass region, the initial decay from the 
plateau ye/ (1 — A) is relevant. To describe this, Eq. (fT6|) 
can be solved by choosing u — and v = 2. The re- 
sulting natural time scale for the initial decay is given by 
h = \fs/ \i\- The long time series has to be matched to 
this time scale. With the same argument as for the liquid 
region, this leads to 



a x (e/e-y -> 00) = a\ (t ) c \e\ 2 \-yt \~ 



(28) 



where a^ is another constant. 



4 The yield process 

While in the previous section, the dynamics of the correla- 
tor around the plateau f c was studied, now its subsequent 
decay is studied. This dynamics is called the a-relaxation 
process in fluids [13j . or the 'yield' process in glasses under 
shear. In this section we derive equations of motion for the 
shear induced final decay (yielding) of <P (t) to zero. First 
we rewrite Eq. (|22j) by introducing a natural time scale: 



Q® 



-t 



£ <M7iT(i7ir, r n (t) 



-cn+1 



n>l 



(29) 



(30) 



3.4.1 Liquid 

In the liquid region, the crossover from the von Schweidlcr 
dynamics to the shear-induced long time dynamics is of 
importance. To obtain the appropriate series expansion 



Now we consider Eq. ([T]). By rescaling time and elim- 
inating the short time dynamics we obtain: 

1 pi 

3(t)-m {£) + —- rh (i - i 1 ) 3 (i') di' = 0, (31) 

dt Jn 
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Fig. 3. Numerically determined master functions of the 'time- 
shear-superposition principle' (dotted lines) 3o (t) (0) and 
$1 (t) (1). The solid lines show numerically obtained transient 
density correlators for e = and 7 = 10~ 7 (red, a), 7 = 1CF 9 
(green, b) and 7 = 1CP 12 (blue, c), plotted as functions of the 
rescaled time t. The plots demonstrate that the rescaled cor- 
relators converge to $q (t) for 7 — > 0, the blue curve (c) is 
already quite close to the master curve (0). The dashed line 
shows <Po (t) + ax |7| c ^i (t) for 7 = 10~ and the same nu- 
merical value for a\ as in Fig. [2] This first order expansion 
already describes quite well the shear-induced decay of the red 
curve (a). The inset demonstrates that the master function 
<&Q (t) (dotted line) can be well approximated by the exponen- 
tial function given by Eq. (|42p (solid line). The curves overlap 
completely. 



*(9 = **® + ,*f ( * ) , (32) 



1+(7*T 




(33) 



Note that 7 does not dependent on 7. Because 7 is a con- 
stant of the order one, we can neglect the influence of e: 



(vi,v 2 ) = {Vi,V$). 



(34) 



With this, the explicit dependence on e and 7 is elim- 
inated. Eq. (|29|) motivates the ansatz of a generalized 
'time-shear-supcrposition principle': 



*(e,7,*)=EK(e,7)r*n(t) 



(35) 



By substituting this in Eq. (|31[) we obtain a recursively 
defined sequence of integro-diffcrcntial equations: 



#„ (t) = J2 $n-n> (t) $n> (t) • 



(38) 



For short rescaled times we require $ (t — ► 0J = f c + 
(1 - f c ) 2 g (f) with Q (t) given by Eq. $25$. Hence we pos- 
tulate: 



#n 



X(e,7) =ai(e,7)|7| c , 
<^o (i -> 0) = f c - (1 - f c ) 2 1, 
(i^0)=-(l-fc) 2 a:^\^)- cn (t) 



— cn+1 



(39) 

(40) 
(41) 



Now we have to show that the postulated short time 
asymptotes represent solutions of theequations of motions 
for t — ► 0. For this we use the ansatz $0 (t — ► 0) = f c +ctoi, 

$n>i (t — ► 0) = a n (i) and neglect (jt) . First we 

,~ 2 
consider n = and neglect all terms of the order [tj . 

The resulting equation is satisfied for arbitrary ao- Now 

we consider n > 1. We can write W n (i) = 2<j> (t) 3 n (i) + 

Y^7=i &n-n' (t) $n' (t) where the first term is of the or- 
der (t) . The sum vanishes for n = 1 and is of the 
order m for n > 2, a sub-leading term. Hence we 
can write: <^ (t — ► 0) = 2& (*) ^n (*) • Similarly, we can 

?) £„< (£') di' is of the 



Itln. — r 



t- 



show that I /„ YZ-=i ■ 

order [tj . We substitute these results in the equa- 

tion of motion and neglect all sub-leading terms. Then the 
resulting equation is satisfied for arbitrary constants a n . 

We remark that integrals like /„ (t') dt' do not ex- 

ist for n > 2/c, but we can formally define objects like 



/o ( ? ) 



' cn+1 di> 



\im S - 



di JS 



Is 



.dt' 



(*") 



-cn+1 



for 



arbitrary n. Such definitions were also implicitly used for 
the series expansions for the glass stability analysis. 

With this we can state: We have derived a closed set 
of equations for <P n (tj which, combined with the appro- 
priate initial conditions, can be solved recursively. The 
solutions are not dependent on e and 7. While the form of 
the correlators is model-dependent, their power-law initial 
decays follows from the universal stability equation ([T0[) . 
It will be these power-law variations that determine the 
exponents in the dependence of the flow curves on shear- 
rate. We remark that n = reproduces the result shown 
in [21] and for e = Eq. ([33)) reduces to the series expan- 
sion presented in [27| . Fig. [3] shows numerical results for 
the master functions $0 (t) and $1 (tj which are consis- 
tent with our results for the glass stability analysis. We 
also observe that the master function <l>o (t) can be well 
approximated by an exponential function: 



^0 (£) = fc exp (-i/io) ■ 



(42) 



d /"' ™ 
$n (t) -fh n {t) + - V m„_„- (i - i') 3 n > (i') di' = 0, 

MJo n'=0 

(36) 

vfP n (i)+v^ n (i) 
m n [t) = -—2 , (37) 

1 + (it) 



The time scale t Q = f c J (1 - f c ) follows from Eq. (j40|) . 

For the discussion of the flow curves we are especially 
interested in n < 2/c where for some finite 5 > the 

integral L <P n (i) di exists. If we assume that J>„ (i) is 
bounded for all i > <5 and all < n < 2/c, then the long 
time limit of rh n (i) is obviously zero. Then by analyzing 
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the Laplace-transforms of the equations motion we can 
easily show: 

Jim 3 n (t) = 0. (43) 

t — *oo 

In addition, we assume for the following that the decay to 
zero is such that <P n (t) is integrable for all < n < 2/c. 

Then for instance, the integral j L <P^ (P) di defines the 
critical yield stress at the glass transition point. 



5 Flow curves: asymptotic regime 

5.1 The yl-formula 

The results of the previous sections, especially Eq. (|35[) 
and Eq. (|39p . allow us to derive an analytical expression 
summarizing many aspects of the flow curves. 

For sufficiently small e and I7I, because of the stretched 
dynamics, the dominant contribution to the integral in 
Eq. ((5]) is determined by the final decay of <P 2 (t) from 
the plateau to zero. Eq. (f35|) and Eq. ([39]) represent an 
asymptotic expansion for the shear-induced final decay of 
<P(t). As demonstrated in Fig. [3] the first few low order 
terms are sufficient to describe the long time asymptote, 
provided that separation parameter and shear rate are 
sufficiently small. We notice that because of Eq. (j40j) and 
Eq. (gT]) the integrals 



m I Y] ®n-n> (t) $n< (*) di 
Jo n >=0 



(44) 



exist for n < [2/c\ = 3 (here [z\ _ denotes the maximum 
integer number n with n < z). Because of the rcscaling of 
time in Eq. ((30|) . the c* are then constants. By taking 
only the integrable low order terms into account, we can 
write down an asymptotic formula for the shear stress: 



<J = ^Z a *n («1 \l\ C ) 7 



(45) 



n=0 



Neglected higher order terms O ( |ai| c |-y| J are subdom- 

inant. Now we use the two-parameter scaling law for a\ 

1 
with the choice J? = lytol 1+a an d define 



A(x) 



ffli {x, 1) 
01 (0,1) 



for some arbitrary argument x. With the definitions a^ = 
(»i(0, !))"<, 

(47) 



1 



c 
2' 



(48) 
(49) 



o 




log 10 (|x|) 

Fig. 4. Numerically obtained data for A (x) (crosses for x < 0, 
plus signs for x > 0) and the fitted functions vl< (x) for small 
negative x (red, 1), AZ, (x) for large negative x (orange, 2), 
/1< (x) for small positive x (green, 3) and At, (x) for large pos- 
itive x (blue, 4), see Appendix lAl for details. With these func- 
tions A (x) can be approximated for arbitrary arguments with 
sufficiently high accuracy. The dotted lines show the leading 
asymptotes with the exponents m" and m! . 



we obtain the yl-formula: 



c(e,7) 



E 

n=0 



A - It* 






(50) 



l7*ol 



(51) 

This is the central result of our work. Eq. (|50|) generalizes 
the results presented in [21] and [57] and summarizes many 
aspects of the flow curves. 

Eq. (fB"0"]) was obtained by integrating the correlators 
obeying the generalized 'time-shear-superposition princi- 
ple'. Their amplitudes depend sensitively on shear rate 
and separation from the glass transition, as could be de- 
termined from the glass stability analysis. One function 
A (e/e-y) = a 1 (e/e-y, 1) ja\ (0, 1) remained capturing this 
sensitive dependence, and in turn now determines the flow 
curves. In the following we discuss the properties of the 
yl-formula in detail. General requirements for the range of 
validity are |e| <?C 1 and |'yto| ^ 1- The <r° are constants, 
which for our model numerically are obtained from fits. 



(46) 5.2 The yl-function 



The central aspect of our work was to analyze the shape 
of yl (x) which we can determine numerically from Q (t) . 
This scalar function of one real parameter determines the 
shape of the flow curves within the range of validity of 
the yl-formula. For instance, the asymptote for x — ► —00 
describes the flow curves in the liquid region while the 
asymptote for x — * 00 describes the scaling of the yield 
stress for e > 0. We can determine three limiting behav- 
iors. 
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-10 
log 10 (y) 

Fig. 5. Overview of the numerically obtained flow curves 
(solid lines) and the yl-formula evaluated numerically (dashed 
lines). The liquid curves in the upper panel are shown for e = 
-10" 7 (red, 1), e = -10 -6 (orange, 2), e = -1(T 5 (green, 3) 
and s = — 10~ 4 (blue, 4). The lower panel shows the glassy 
curves for e = (red, 1), e = 10 -5 (orange, 2), e = 10 -4 
(green, 3) and e = 10 -3 (blue, 4). Crosses mark the points with 
|e| = e-y. The natural upper boundary for the shear rate, where 
the range of validity of the A- formula is left, is also indicated, 
see Eq. (|59p . For e < 0, the natural lower limits for the shear 
rates described by Eq. (|62|l . below which the yl-formula does 
not describe the flow curves, are marked by diamonds. 



By using Eq. fl26 
A(x- 



Eq. (HZJ) and Eq. J28J) we identify 
-oo) = -/I" \x\ m " , (52) 



(53) 

(54) 
(55) 



"~ l-b' 1 
A(0) = 1, 

/l(x^oo)=/l+x m '. 

We determine A^ and /1+ numerically. Fig. |4] shows an 
overview of the properties of the /1-function. A more de- 
tailed discussion can be found in Appendix [AJ 

5.3 Discussion of the flow curves 

Fig. \E\ shows an overview of the numerically obtained flow 
curves and the corresponding numerical results for the yl- 
formula. While the glass flow curves exhibit an upward 
curvature only, the fluid curves show a characteristic S- 
shape, where the initial downward curvature changes to 
an upward one for increasing shear rate. For positive sep- 
aration parameters the range of validity of the yl-formula 
is given by |e| <C 1 and \jto\ <C 1. These two requirements 
ensure that Q (t) describes the dynamics of $ (t) with a 
sufficiently high accuracy. For sufficiently small negative 
separation parameters the yl-formula is valid in finite shear 
rate windows only, as it does not reproduce the linear 
asymptotes for low shear rates. We find a precise criterion 
for the range of validity later when we discuss the limiting 
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Fig. 6. Numerically obtained curve for the yield stress (solid 
line). The dashed lines show the corresponding power series 
given by Eq. (|56p evaluated to first order (red, 1), second or- 
der (green, 2) and third order (blue, 3). The natural upper 
boundary for the range of validity of Eq. (|56|l , respectively Eq. 
(1571). is also indicated. 



In the framework of the yl-formula, it is easy to under- 
stand the mathematical structure of the singular behavior 
of the yield stress at e = 0: We choose some fixed e = ± |e| 
and consider the limit 7 — » 0. Then for x = s/s-y we obtain 
x — > ±00. The singular behavior is then due to the differ- 
ent asymptotes of A (x) for x — > ±00 while yl (x) itself is 
a smooth function for all x. 

While Fig. [5] presents an overview of the flow curves, 
detailed discussions of the various limits are presented 
next. 



5.3.1 Yield stress 



By using Eq. (|55|) we obtain an asymptotic formula for 
the yield stress which reproduces a result in [2"4"j : 



a + (e) 



E 

71 = 



r +~m n 



(56) 



The constants are defined by ct+ = cr" (yl+ )". We deter- 
mine CTq = a + [e = 0) by extrapolating the flow curve for 
£ = 0. Then we use a non-linear curve fitting to determine 
°n>i- The quantity ctq is the dynamic yield stress at the 
transition, also denoted by er+; numerically we verify that 
it agrees with the integral c+ = 7 L $Q(t)di, as implied 
by the time-shear superposition principle Eq. (|35[) . Equa- 
tion (f56|) describes the yield stress for sufficiently small 
positive separation parameters, see Fig. [HI and can be 
rewritten as 



ct+ o) = 0-0" Yl c « ( e / e *) m 



(57) 



n=0 



with Cq = 1. With the choice c| = 1 we can identify 

1 
£ * = ( (J o/ a i)~™ T as a natural scale for the separation 
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Fig. 7. Numerically obtained flow curve for s = (solid line) . 
The dashed lines show the corresponding power series given by 
Eq. (|58[) evaluated to first order (red, f), second order (green, 
2) and third order (blue, 3). The natural upper boundary for 
the range of validity of Eq. ([58}, respectively Eq. (|59|) . is also 
indicated. 



parameter in our asymptotic expansion. To keep inside the 
range of validity of the A- formula we have to require |e| -C 
e*. For the constants we obtain c^ >2 = (&n / a o) e ™ "• 



5.3.2 Critical curve 

By using Eq. (|54| we obtain an asymptotic formula for 
the critical flow curve reproducing the results shown in 
[24] and [27], the generalized Herschel-Bulklcy law: 



L£J_ 

<r( e = 0, 7 ) = Yl ^l^ol 



(58) 



We determine the constants directly from the fitted values 
for cr+ by using er+ = a° (/1+ )". We observe that for e = 
Eq. ([58]) describes the flow curve correctly for sufficiently 
small shear rates, see Fig. [7] Eq. ([55]) can be rewritten as 



(59) 



a (e = 0, 7) = o-+ J^ c n I7/7* 



With the choice Cq = C\ = 1 we can identify <t+ = erf, as 

1 j_ 

the critical dynamic yield stress and 7* = (to) ( <T o/ <7 i) m 
defines a natural scale for the shear rates in our asymp- 
totic expansion [8]. To keep inside the range of validity 
of the A- formula we have to require |7| ^C 7*. For the 
constants we obtain c„>2 = (c^/cg) ( A f*to) mn - 



5.3.3 Liquid 



Eq. ([52]) leads us to a formula for negative separation pa- 
rameters: 

cr(£<-£^,7)= ^ (-l)"cr 



— 1 iffl n 



7*0 1 



(60) 




-20 -19 

log 10 (Y) 

Fig. 8. Numerically obtained flow curve for e — — 10 -9 (solid 
line). The dashed lines show the corresponding power series 
given by Eq. (|60p evaluated to first order (red, 1), second order 
(green, 2) and third order (blue, 3). The shear rate with e = 
— e^ is marked by a cross. The natural lower limit A f m in for 
the shear rate defined by Eq. (|62p is marked by a diamond. A 
fitted regular Taylor expansion to the order 7 s for the region 
below 7 m in is also included (dotted line). We remark that we 
need extremely small e and 7 to verify Eq. (|60[) numerically. 



7 



(61) 



We determine the constants directly from the values for 
Og by evaluating a~ = cr° (-A~) n . We observe that for 
sufficiently small separation parameters Eq. ([60]) describes 
the flow curves correctly in finite shear rate windows where 
the flow curves are located close below the critical yield 
stress value, see Fig. [8] 

A necessary criterion for the validity of Eq. (|60p is 

b 1 

that the matching time t& cx \e\ 1_6 |7<o| 1 ~ b lies inside 
the range of validity of the nonlinear stability equation, 
see Eq. ([TO]) . This criterion is equivalent to the statement 
that Q {t) converges to the long time asymptote —t/r^ un- 
der the condition \Q (t)\ -^ 1. This consideration leads to 
lower limits for the shear rates, below which the yl-formula 
does not describe the flow curves. They are marked by di- 
amonds in Fig. [5] and Fig. [8] The requirement r;, <^i re- 
loads to I7I > 7 mm with 



(f t 



- 1 6) 



(62) 



Below this limit, a regular Taylor expansion holds which 
is also shown in Fig. [8] 



5.3.4 Yielding glass 

Another interesting aspect is that for e > we observe 
numerically a linear regime for small shear rates: 



a(e > 0,7^0) = a + + a x I7I 



(63) 



This result enables us to identify the next-to-leading order 
asymptote of the A- function for large positive arguments, 
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Numerically obtained flow curves from which the yield 



stress is subtracted (solid lines) for e = 10 
(orange, 2), e = 10~ 5 (green, 3) and e — 10 



7 (red, 1),£ = 10" 6 
" 4 (blue, 4). Crosses 
mark the points with e = Ej. The dotted lines show the fitted 
linear functions ai |7|. The inset demonstrates that the fitted 
coefficients (crosses) obey ai tx e^^ (solid line) for sufficiently 
small e. 

see Appendix [A] With this we can predict a scaling law 
for the coefficient <j\ for small e: 



<7i oc £ 



c= 



1 — m 



(64) 
(65) 



These formulae are consistent with our numerical results 
shown in Fig. [9] Note that Eq. ([63]) defines a function 
which is not analytic in 7. 



5.4 The pseudo power law 

For e < we observe that log 10 (cr) as a function of 
log 10 (7) shows an inflection point defined by 



d 2 (login (a)) 
rf(!ogio(7)) 2 



0. 



(66) 



Obviously, this statement is also true for any (toy-) model 
which shows shear-thinning. But then in some finite shear 
rate windows the flow curves can be approximated by the 
corresponding inflection tangents. The slopes p of the in- 
flection tangents can be interpreted as exponents occur- 
ring in some pseudo power laws: 



a oc 7 P . 



(67) 



Fig. [TU] shows some examples. We observe that the numeri- 
cally evaluated /1-formula also describes the neighborhood 
of the inflection point correctly for sufficiently small e. As 
shown in the plot, in this region — e~ t < e < is satisfied, 



hence the A- function can be described by Eq. (|76)) . see 
Appendix [A] for details. By substituting Eq. (|76[) in Eq. 
(|5D|) we see that the resulting formula does not represent 
a real power law of the type of Eq. (|6T|) . 



iog 10 (y) 



Fig. 10. The upper panel shows numerically obtained flow 
curves (solid lines) for e = — 10 -3 (red, 1), £ = — 10~ 2 (green, 
2) and e = — 10 _1 (blue, 3). The dotted lines show the cor- 
responding inflection tangents, with exponents p — 0.16, 0.35, 
and 0.63 from left to right. The dashed line shows the nu- 
merically evaluated A- formula for e = — 10~ 3 . The shear rate 
with e — —e-y is marked by a cross. The lower panel shows the 
corresponding results for the viscosity. 



We conclude that, in the framework of our asymp- 
totic expansion, there is no real exponent p. Eq. (|67|) is 
a trivial artifact of the double logarithmic plot of the 
flow curves which then necessarily show inflection points 
for any model describing shear thinning. Rather, the flow 
curves on the fluid side exhibit a characteristic S-shapc. 
While this shape is rather apparent when plotting stress as 
function of shear rate, plotting the same data as viscosity 
as function of shear rate hides it. The vertical axis gets 
appreciably stretched, and ruling out the pseudo power 
law becomes more difficult using experimental data. 



6 Flow curves: intermediate regime 

The correlator <P(i) from Eq. (|31[) determines the struc- 
tural contribution to the stationary stress, where micro- 
scopic, initial-dccay-ratc--r-depcndcnt contributions are ne- 
glected. It contains the complete structural dynamics on 
time scales longer than in, and starts out with the critical 
power law according to Eqs. ([9]) and (fTTj) . Expansion (|35|) 



captures its behavior for rescaled times t 3> Je~^ only, be 



cause the power series (P^|) . on which expansion ([55)1 is 
based, at this scale crosses over to expansions on shorter 
time scales; the matchings are discussed in Sect. 3.4. The 
divergence of the expansion coefficients # n (i) for i — > 
thus is only apparent, and becomes cut off at t = 0{*fe^). 
Let us define a correlator <Z>(t), where the (dominant) 
divergent initial variations arc subtracted: 

*(t) = *(t) - (l - fcf V^ G(i/V^) °( l - *) . ( 68 ) 

where Q is the series from Eq. (P^|) valid for Je^ <C 
i « 1, only. The new correlator possesses an expansion 
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Fig. 11. Critical flow curves for three different transitions 
and exponent parameters A rescaled using the critical yield 
stress a^ and the scale 7, to agree with the Herschel-Bulkley 
law from Eq. (J7TJl with exponents m = 0.200 (A = 0.976, red, 
label 1), m = 0.489 (A = 0.707, green, label 2), and m = 0.542 
(A = 0.577, blue, label 3); curves 2 and 3 are shifted upwards by 
I and 1, respectively. Solid lines are the numerical flow curves, 
dotted lines the generalized Herschel-Bulkley law from Eq. (|59[) 
which holds for 7 < 7*, and dashed lines give Eq. (|71|l which 
holds for 7 ~ 7*. The inset shows the multiplied derivatives 
7 ~ m da/dj, where the corresponding scales 7* are marked. 
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Fig. 12. Reduced flow curves for a core-shell dispersion at 
an effective volume fraction of (fieff — 0.580. Here Rh denotes 
the hydrodynamic radius and Do the self diffusion coefficient of 
the colloidal particles; k B T is the thermal energy. The solid line 
(red) shows the result for the fitted F^ -model with «| = 2.0. 
The fitted parameters are: e = —0.00042, 7 C = 0.14, v a — 
70k B T/R 3 H , r = 80Do/R 2 H and i loo = 0.3Q4k B T/R H D . The 
dashed line shows the corresponding result for the Tl-formula. 
The dotted line shows the inflection tangent of the numerically 
determined flow curve with a slope of p = 0.12. The inset shows 
the corresponding results for the viscosity. 



like Eq. ([55)) with coefficient functions <£„(£) that behave 
more regularly at short rescaled times. It thus can be inte- 
grated over rescaled time t. The Heaviside-step function in 
Eq. ([68]) serves as remembrance that the series ([22]) only 
describes the initial variation of the yield scaling functions 
S. The sum, 

#(£) = #(t) + (1 - f c f Je~ G(i/^e-) 0(1 - i) , (69) 



where now the complete /3-correlator Q from Eq. ([10]) is 
taken, thus describes the solution of Eq. ([31| down to its 
limit of validity, |'yio I ^ t- 

In the limit of small shear rates, |'y£o| ^ 1) the stress 
from Eq. ([5]) using this rewriting of $ becomes: 



7 



(*(£) + (1 - f c f 



7n\2 



G(i/^e-)0(l-i)) di 



7 / <P 2 (i)di +h, 



(7 £7 



1/ 



3/2n 



G{t)dt + 0« z ) (70) 



where t = t/ \f£^\ the constant abbreviates h a = 27 f c (1 — 

Icf- 

The result just obtained agrees with Eq. ([59)) in the 
asymptotic limit 7 — ► 0. This holds because the scaling 
function Q(t) diverges for long rescaled times i, and thus 
the behavior of the second integral in Eq. ([70)) at the up- 
per limit, which shifts to infinity in this case, gives the 
leading contributions. Yet, the scaling function Q remains 
close to zero for rescaled times t, which correspond to the 
region where the full correlator stays close to the plateau 
/ c ; see Fig.Q] For not too small shear rates, which we es- 
timate based on Eq. ([5"9")) by |7|^7*, the second integral 



in Eq. (j70[) thus does not depend on 7, and we derive a 
Herschel-Bulkley law for intermediate shear rates: 

a(% « | 7 | « r,e = 0) = a {) + a x \^t f + O(| 7 t | 3 " V2 ) 

(71) 
with constant amplitudes 



'0 



i/v^T 



Q{i) di 



The constant <7o is not the actual yield stress, cr + , which is 
obtained in the limit of vanishing shear rate, a + = a(j — > 
0), but is larger, because <£>o(£) remains longer close to the 
plateau value f c than <Pq (t) , which enters in the definition 
of <t + ; this follows because <Po(t) lacks the initial term 
linear in t. 

Fig. QT] verifies the Herschel-Bulkley law for three dif- 
ferent exponents m = jt^, obtained by choosing three 
different glass transition points characterized by different 
exponent parameters A introduced in Sect. 2.1. The expo- 
nents m vary by almost a factor of 3, and the window in 
shear rates, where Eq. 1(71)) holds, shifts appreciably with 
A. The plot of j 1 ~ rn da/d , j versus shear rate shows an (al- 
most) constant plateau in the 7-window, where Eq. (|71|) 
holds. 



7 Comparison with experiments 

Fig. [T2] and Fig. Q2] show experimental data recently ob- 
tained by Siebenbiirger et al. [35] on polydisperse dis- 
persions of thermosensitive core-shell particles [10) . In all 
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Fig. 13. Reduced flow curves for a core-shell dispersion at 
an effective volume fraction of </> e // = 0.629; quantities as de- 
fined in the caption of Fig. 1121 The solid line (red) shows the 
result for the fitted Ff^'-model with u| = 2.0. The fitted pa- 
rameters are: e = 0.000021, j c = 0.16, v a = 115k B T/R :i H , 
r = 120Do/R 2 H and T] x = 0A31k B T/R H D . The dashed line 
shows the corresponding result for the /1-formula. The dotted 
line shows the fitted Herschel-Bulkley law given by Eq. (|71[1 
with the analytically calculated exponent rh = 0.489. 



cases stationary states were achieved after shearing long 
enough, proving that ageing could be neglected even for 

glassy states. Along the lines of the work [8], the F^ - 
model was used to fit both the flow curves and the lin- 
ear viscoelastic moduli simultaneously, which strongly re- 
stricted the choice of fit parameters. Here we only show 
two representative results for the flow curves, comparing 
the asymptotic results. 

Fig. [12] shows the result for a liquid-like flow curve 
where the /1-formula holds for approximately four decades. 
The pseudo power law resulting from the inflection tan- 
gent of the flow curve holds for approximately two decades 
within the range of validity of the /1-formula. 

Fig. fT3l shows the result for a flow curve, where a small 
positive separation parameter was necessary to fit the flow 
curve and the linear viscoelastic moduli simultaneously. 
The data fall below the fit curves for very small shear 
rates, which indicates the existence of an additional de- 
cay mechanism neglected in the present approach [8|35j . 
Again, the /1-formula describes the experimental data cor- 
rectly for approximately four decades. For higher shear 
rates, the Herschel-Bulkley law given by Eq. ([71"]) can be 
fitted in a window of approximately two decades. We have 
derived Eq. (|7Tj) for e = 0, nevertheless the fit is possible 
for two reasons: In the present case, the separation pa- 
rameter is sufficiently small, and for higher shear rates 
the e-dependence of the flow curves becomes negligible. 



8 Conclusion 

In the present work we have analyzed the shear flow be- 
havior of dense colloidal dispersions in the framework of 



the schematic F^ -model based on the mode-coupling 
theory. This model describes the transition from a shear- 
thinning liquid to a yielding glass [24] . While in the liquid 
state the steady state shear stress decays linearly to zero 
for small shear rates, in the glassy state the shear stress 
arrests on a finite dynamic yield stress value in the zero- 
shear limit. In the liquid thus a linear response regime 
characterized by a Newtonian viscosity exists. The viscos- 
ity diverges at the glass transition according to Eq. ([20]) 
[13j , and carries a regular shear rate dependent correction 
of order j 2 . The glass, on the other hand, never exhibits a 
regular linear response in shear rate. For high shear rates 
the shear stress depends linearly on the shear rate for both 
cases. 

Close to the glass transition and for small shear rates, 
we studied the model using asymptotic expansions. The 
derivation of a quite general asymptotic formula summa- 
rizing many aspects of the flow curves, viz. the steady 
state shear stress as a function of the shear rate, was the 
central result of our work. This formula provides a deeper 
understanding of the shapes of the flow curves. The scal- 
ing of the yield stress and its singular behavior at the 
glass transition point can easily be obtained by studying 
special asymptotic limits. The yield stress a + jumps dis- 
continuously from zero in the fluid to a finite value at 
the glass transition, and grows with the power law series 
from Eq. ([57)1 deeper into the glass. For all states (deep) 
in the glass, the stress increases (non-analytically) with 
the absolute value of the shear rate, see Eq. ([63]) . At the 
glass transition point our formula reduces to the general- 
ized Herschel-Bulkley law Eq. ([59]) known from previous 
studies [24j and [27] . We have also demonstrated in Fig. \W\ 
that the apparent power law behavior of the flow curves, or 
of the viscosity, in the liquid state for intermediate shear 
rates is only an artifact of the inflection point occurring 
in a double logarithmic plot. This inflection point neces- 
sarily occurs for any model including shear thinning. The 
Herschel-Bulkley law of Eq. ([7T]) . which holds in a win- 
dow of intermediate and not too small shear rates, should 
be quite accessible experimentally Importantly, the expo- 
nent rh is connected to other exponents characterizing the 
material properties, and can thus be inferred from other 
measurements. Also, it can be calculated using Eq. ([T2]) 
from first principles starting with the particle interactions 
and using mode coupling theory [13] . 

In previous studies, the F^ -model was successfully 
applied to describe experimental data for both the steady 
state shear stress and the linear viscoelasticity, see [8] and 
[TO] . As shown in [27] and [5] , the model also describes sim- 
ulation data for the steady state shear stress with a high 
accuracy. Recently, the Herschel-Bulkley law of Eq. ([7T]) 
with the predicted exponent m HS = 0.473 was observed in 
experiments on dispersions of colloidal hard spheres |29j . 
However, only the recent experiments of Siebenbiirger et 
al. [33] provide access to shear rates and distances from 
the glass transition points which arc sufficiently small for a 
test of our asymptotic description of the flow curves. Our 
asymptotic formulae describe these experiments success- 
fully in shear rate windows of several decades. Especially, 
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Fig. 14. The next-to- leading order asymptotes of A (x). The 
crosses are numerically obtained data for A (x). The two upper 
panels show linear fits for small \x\. The left lower panel shows 
the power law for x —> — oo with the fitted exponent £. The 
right lower panel shows for x — > oo the fit with the analytically 
calculated exponent £. 



the pseudo power law is clearly observable in the experi- 
mental data. 
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The constant \ is defined by A (— x) = 0. We obtain \ ~ 2. 
We choose N = 3 and determine A^ and A~ by non-linear 
curve fittings. For positive separation parameters, we have 
identified 

A (x -► 0+) = 1 + A+x, (77) 



A{x 



oo) = A^x™' 



A[a 



-C 



In addition, by using Eq. (|63|) we have identified 

1 — m 



( 



(78) 



(79) 



Inspired by these results, we have found that A (x) can be 
well approximated by 



A(x) 



A+(x) 




At (x) = l + 5>+x". 



(80) 



(81) 



(82) 



n=l 



Again, we choose N = 3 and determine A^ and A+ by 
non-linear curve fittings. The fact that we obtain similar 
fitted values for A£ and A~ is an indication for a smooth 
behavior of A (x) at x = 0. With the presented formulae 
we are able to describe A (x) for arbitrary arguments. Fig. 
@]and Fig. [Til summarize the results for the A- function. 



B Numerical parameters 



A Discussion of the /l-function 

In this section we present some numerical results for the 
next-to-leading order asymptotes of A (x) . For negative 
separation parameters, we have identified 



A (x^ -*>) = - A^\x\ m "+ Af\x\- 



A (x -> (T 



1 + Ajx, 



(72) 
(73) 



with an empirical exponent £. Inspired by this, we have 
found that A (x) can be well approximated by using the 
ansatz 



A (x < 0) 



A > (x) , |x| > x 
A~(x) ,0<|x|< x 



A ! 



A-(x) = -A^\x\ m +J2^i 



-«» 



/V 



A-(x)=l + J2 A nX n - 



(74) 

(75) 
(76) 



The numerical values for the parameters, which are rele- 
vant for our work, are summarized in the Tables [TJ|4] 



vl 


0.828 


A 


0.707 


to 


0.426 


T„ 


0.240 


V2 


2.00 


C W 


0.586 


To 


1.71 




0.0858 


fc 


0.293 


7 


0.595 


to 


0.586 







(-V) 



Table 1. The parameters for the schematic F^ -model. We 
have chosen v^ = 2.00. The time scales to, ^o and f v were cal- 
culated numerically. All other quantities were calculated ana- 
lytically. 



ra=l 
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Table 2. The numerical values for the exponents. The expo- 
nents £ and £ were first identified numerically. In addition, an 
analytical expression could be verified numerically for £ which 
was used to calculate it. All other exponents could be derived 
analytically. Hence we can state that these exponents and £ 
are only dependent on the exponent parameter A, in this sense 
they are universal. This statement could not be clearly verified 
for^. 
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Table 3. Numerically determined amplitudes for the lead- 
ing asymptotes and the empirical correction terms for the A- 
function. The fact that we obtain similar fitted values for A„ 
and AX is an indication for a smooth behavior of A (x) at x = 0. 
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Table 4. Numerically determined amplitudes for the flow 
curves. The only fitted parameters are a„ from which the other 
quantities were calculated. 
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